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The cell cytoskeleton is a striking example of "active" medium driven out-of-equilibrium by ATP 
hydrolysis 1 . Such activity has been shown recently to have a spectacular impact on the mechan- 
ical and rheological properties of the cellular mediu m 1 2 * 3 * 4 * 5 * 6 * 7 * 8 * 9 * 10 !, as we ll as on its transport 
propertie d 11 * 12 * 13 * 14 !: a generic tracer particle freely diffuses as in a standard equilibrium medium, 
but also intermittently binds with random interaction times to motor proteins, which perform ac- 
tive ballistic excursions along cytoskeletal filaments. Here, we propose for the first time an analytical 
model of transport limited reactions in active media, and show quantitatively how active transport can 
enhance reactivity for large enough tracers like vesicles. We derive analytically the average interaction 
time with motor proteins which optimizes the reaction rate, and reveal remarkable universal features 
of the optimal configuration. We discuss why active transport may be beneficial in various biological 
examples: cell cytoskeleton, membranes and lamellipodia, and tubular structures like axons^. 

Various motor proteins such as kinesins or myosins are able to convert the chemical fuel provided by ATP into 
mechanical work by interacting with the semiflexible oriented filaments (mainly F-actin and microtubules) of the 
cytoskeleton^. As many molecules or larger cellular organelles like vesicles, lysosomes or mitochondria, hereafter 
referred to as tracer particles, can randomly bind and unbind to motors, the overall transport of a tracer in the cell 
can be described as alternating phases of standard diffusive transport, and phases of active directed transport powered 
by motor proteins' 1 * 15 * 16 * Active transport in cells has been extensively studied both experimentally, for instance by 
single particle tracking methods' 11 * 12 !, and theoretically by evaluating the mean displacement of a tracer-^^, or 
stationary concentration profiles^. 

On the other hand, most of cell functions are regulated by coordinated chemical reactions which involve low 
concentrations of reactants (such as ribosomes or vesicles carrying targeted proteins) , and which are therefore limited 
by transport. However, up to now a general quantitative analysis of the impact of active transport on reaction kinetics 
in cells, and more generally in generic active media, is still missing, even if a few specific examples have been tackled^. 
In this letter, we propose an analytical model which allows us to determine for the first time the kinetic constant of 
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transport limited reactions in active media. 

The model relies on the idea of intermittent dynamics introduced in the context of search 
processe d 19 * 20 * 21 * 22 * 23 * 24 * 25 * 26 * 27 ! We consider a tracer particle evolving in a d-dimensional space (in practice 
d = 1,2,3) which performs thermal diffusion phases of diffusion coefficient D (denoted phases 1), randomly 
interrupted by ballistic excursions bound to motors (referred to as phases 2) of constant velocity v and direction 
pointing in the solid angle a-> v (see figure la). The distribution of the filaments orientation is denoted by p(ui v ), and 
will be taken as either disordered or polarized (see figure la), which schematically reproduces the different states of 
the cytoskeletoiW The random duration of each phase i is assumed to be exponentially distributed with mean Tj. 
The tracer T can react with reactants R (supposed immobile) during free diffusion phases 1 only, as T is assumed to 
be inactive when bound to motors, which is realized for instance when reactants are membrane proteins (see figure 
lb,c). Reaction occurs with a finite probability per unit of time k when the tracer-reactant distance is smaller than 
a given reaction radius a. In what follows we explicitly determine the kinetic constant K of the reaction T + R — > R. 

We now present the basic equations in the case of a reactant centered in a spherical domain of radius b with 
reflecting boundary. This geometry both mimics the relevant situation of a single target and provides a mean field 
approximation of the general case of randomly located reactants with concentration c = a d /b d , where b is the typical 
distance between reactants. We start from a mean field approximation of the first order reaction constant! 2 ^ and write 
K = l/(£), where the key quantity of our approach is the reaction time (t) which is defined as the mean first passage 
time (MFPT) 29 * 30 * of the tracer at a reactant position uniformly averaged over its initial position. For the active 
intermittent dynamics defined above, the MFPT of the tracer at a reactant position satisfies the following backward 
equatiorP^ 

!DA r ii + — / (i 2 - h)p(uj v )duj v - fcl (r)ii = -1 
v « 
v-V,ta (*a-*i) = -l 
T2 

where t\ is the MFPT starting in phase 1 at position r, and t 2 is the MFPT starting in phase 2 at position r with 
velocity v. I a is the indicator function of the ball of radius a. As these equations ([TJ are of integro-differential type, 
standard methods of resolution are not available for a general distribution p. 

However, in the case of a disordered distribution of filaments (p(uj v ) — l/f2<j, where Qd is the solid angle of the 
d-dimensional sphere) , one can use a generalized version of the decoupling approximation introduced irPSl t obtain 
a very good approximate solution of equations 0, as described in the Methods section. We present here simplified 
expressions of the resulting kinetic constant by taking alternatively the limit k — » oo, which corresponds to the ideal 



case of perfect reaction, and the limit D — > which allows us to isolate the k dependence. 

We first discuss the d — 3 disordered case (see figure la), which provides a general description of the actin cytoskele- 
ton of a cell in non polarized conditions, or of a generic in vitro active solution. An analytical form of the reaction 
rate K 3 d is given in the Methods section, and plotted in figure 2a, b. Strikingly, K 3 d can be maximized (see figure 
2a, b) as soon as the reaction radius exceeds a threshold a c ~ D/v for the following value of the mean interaction time 
with motors: 

r°f d = — * 1.078- (2) 

where Xq is the solution of 2tanh(ai) — 2x + xtanh(x) 2 = 0. The T\ dependence is very weak, but one can roughly 
estimate the optimal value by r"^ ~ 6D/v 2 . This gives in turn the maximal reaction rate 

co %/3(x -tanh(a;o)) 

x% ' (3) 

so that the gain with respect to the reaction rate K^ d in a passive medium is G$d = K^/K^ ~ Cav/D with C ~ 0.26. 

Several comments are in order, (i) First, p ^ d neither depends on D, nor on the reactant concentration. A similar 
analysis for k finite (in the D — > limit) shows that this optimal value does not depend on k either, which proves that 
the optimal mean interaction time with motors is widely independent of the parameters characterizing the diffusion 
phase 1. (ii) Second, the value a c should be discussed. In standard cellular conditions D ranges from ~ 1CP 2 ^m 2 .s _1 
for vesicles to ~ 10 [im 2 .s~ 1 for small proteins, whereas the typical velocity of a motor protein is v ~ 1 /xm.s -1 , value 
which is widely independent of the size of the cargcW. This gives a critical reaction radius a c ranging from ~ 10 nm 
for vesicles, which is smaller than any cellular organelle, to ~ 10 \xm for single molecules, which is comparable to the 
whole cell dimension. Hence, this shows that in such 3-dimensional disordered case, active transport can optimize 
reactivity for sufficiently large tracers like vesicles, as motor mediated motion permits a fast relocation to unexplored 
regions, whereas it is inefficient for standard molecular reaction kinetics, mainly because at the cell scale molecular 
free diffusion is faster than motor mediated motion. This could help justifying that many molecular species in cells are 
transported in vesicles. Interestingly, in standard cellular conditions t% 3^ is of order 0.1 s for a typical reaction radius 
of order 0.1 fim. This value is compatible with experimental observations^, and suggests that cellular transport is 
close to optimum, (iii) Last, the typical gain for a vesicle of reaction radius a > 0.1 /im in standard cellular conditions 
is G^d > 2.5 (see figure 2a, b) and can reach G^d ^ 10 for faster types of molecular motors like myosins (v ~ 4/zm.s -1 , 
se e 1 1 ! 11 * ), independently of the reactant concentration c. As we shall see below the gain will be significantly higher in 
lower dimensional structures such as axons. 



We now come to the d = 2 disordered case (see figure lc) . Striking examples in cells are given by the cytoplasmic 
membrane, which is closely coupled to the network of cortical actin filaments, or the lamellipodium of adhering cells^. 
In many cases the orientation of filaments can be assumed to be random. This problem then exactly maps on the 
search problem studied irf22l, where the reaction time was calculated. This permits to show that as for d = 3, the 
reaction rate K-^d can be optimized in the regime D/v <C a <C b. Remarkably, the optimal interaction time t% \d takes 
one and the same value in the two limits k — > oo and D — > : 

^^(Ml/^-lp) 1 ^, (4) 

which indicates that again r^d does not depend on the parameters of the thermal diffusion phase, neither through 
D nor k. In the limit k — » oo one has t?^, = n wiV^ffi , . and the maximal reaction rate can then be obtained : 

1,2a 2v* 21n(l/vc) — 1 

K™ d ~ ^ (5) 

Comparing this expression to the case of passive transport yields a gain G^d = K^d/K^d — o.vy/]n(l/y/c) / (4D). As in 
the d — 3 case, this proves that active transport enhances reactivity for large enough tracers (with a critical reaction 
radius a c ~ D/v of the same order as in the d = 3 case) such as vesicles. However, here the gain G<id depends on the 
reactant concentration c, and can be more significant : with the same values of D, v and a as given above for a vesicle 
in standard cellular conditions, and for low concentrated reactants (like specific membrane receptors) with a typical 
distance between reactants b > lO/im, the typical gain is G2d ^ 8, and reaches 10 for single reactants (like examples 
of signaling molecules). 

The case of nematic order of the cytoskeletal filaments, which depicts for instance the situation of a polarized celP} 
can be shown to be equivalent in a first approximation to the 1-dimensional case, which is exactly solvable (see figure 
la,b). The d = 1 case is also important on its own in cell biology as many 1-dimensional active structures such as 
axons, dendrites, or stress fibers^ are present in living cells. As an illustration, we take the example of an axon, filled 
with parallel microtubules pointing their plus end in a direction e. We consider a tracer particle interacting with both 
kinesins ("+" end directed motors, of average velocity ve ) and dyneins ("-" end directed motors, of average velocity 
—ve) with the same characteristic interaction time r 2 (see figure lb). For this type of tracer, the MFPT satisfies 
equations |TJ) with an effective nematic distribution of filaments p(oj v ) = §(#(v — e ) + ^( v + e ))- The reaction rate 
Kid is obtained exactly in this case (see Methods), and is maximized in the regime D/v <C a <C b for the following 
values of the characteristic times (see figure 2c,d) 

Opt _ J_ D Opt _ 1 a 

Tl > ld ~ 48 T ^ ld ~ ^3^/2' W 



for k — ► oo. The maximal reaction rate K™ d is then given by 

Kid- 2a > (7) 

and the gain is Gid = K™ d /K^ d ~ au/(2v / 3-Dc 1 / 2 ), which proves that active transport can optimize reactivity as 
in higher dimensions. Very interestingly the c dependence of the gain is much more important than for d — 2,3, 
which shows that the efficiency of active transport is strongly enhanced in 1-dimensional or nematic structures at 
low concentration. Indeed, with the same values of D, v and a as given above for a vesicle in standard cellular 
conditions, and for a typical distance between reactants b > 100/im (like low concentrated axonal receptors), one 
obtains a typical gain G\d > 100 (see figure 2c, d). In the limit of finite reactivity (k finite and D — > 0) one has 
T ° P id = \/c i /'vk(\/12c) 1 / A and the same optimal value (JsJ) of T^ d . As in higher dimensions T^ d depends neither 
on the thermal diffusion coefficient D of phases 1, nor on the association constant k, which shows that the optimal 
interaction time with motors T 2 pt presents remarkable universal features. Furthermore, our approach permits an 
estimate of compatible with observations in standard cellular conditions, which suggests that cellular transport 
could be close to optimum. 

Methods 

The approximation scheme to solve the integro-differential equations ([l| relies on the auxiliary function 

*(r) = ^- J t 2 dw v , (8) 

and on the following decoupling hypothesis: 

v 2 

(viVjta)^ — (viVj) Wv {t2) Uv = — %s(r). (9) 
Similar arguments as provided irP^l then leads to the diffusion-like equation 

DAs(r)-—(s(r)-t 1 ) = -l (10) 

1~2 



where D = v Tijd. After rewriting equation)!!]) as 



DAh + -(s(r) - t x ) - kl a (r)h = -1, (11) 



equations! 10 1 and (111 provide a closed system of linear differential equations for the variables s and t\, whose 



resolution is tedious but standard. For d = 3, we obtain in the limit of perfect reaction k — > oo and low density a -C b: 

3 ai D (£>riT(l + a x a) + D(ti + T 2 )(a 2 a - T)) 
b 3 (T + a 1 a 2 D(T 1 +T 2 )) 
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where ol\ = \/^Dt\, 0*2 = 3/ (1^2) and T = tanh(a!2&). This decoupling assumption has been controlled numerically 
for a wide range of the parameters for d = i-PH, and is shown here to be satisfactory also for d — 3 (see figure 2a, b). For 
d = 1 the decoupling approximation is exact and yields after straightforward calculations an explicit though hardly 
handleable form of the reaction time. In the regime D/v <C a <C b, we obtain in the limit k — > 00 the simple form 




where s — Optimization is then performed using standard methods of functional analysis. 
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Figure legends: 

FIG. 1: Model of reaction kinetics in active media, and examples of low dimensional structures in biological cells, a, 
The reactant alternates thermal diffusion phases (regime 1 in red) of mean duration t\ and diffusion coefficient D, and 
ballistic phases of velocity v powered by molecular motors (regime 2 in blue) of mean duration r 2 . The cytosqueletal 
filaments (in black) are here in a disordered state. The polarized nematic state would correspond to parallel filaments, 
and is equivalent in a first approximation to a 1 -dimensional situation (see b) with the same concentration c\d = 
C3d — a 3d/^3d an d an effective reaction radius aid = a3dCg^ 3 . Molecular motors are not represented, b, Tubular 
structures in cells such as axons and dendrites (d = 1). c, Planar structures such as membranes and lamellipodia 
(d=2). 

FIG. 2: Optimization of the reaction rate. a,b, Gain of reactivity due to active transport G^d for d = 3 as a function 
of T2 for different values of the ratio b/a (logarithmic scale). The analytical form obtained in the Methods section 
(plain lines) is plotted against numerical simulations (symbols) for the following values of the parameters (arbitrary 
units): a = 1.5 (brown), a — 4.5 (red), a — 7.5 (blue), a — 10.5 (green), a = 15 (yellow), with t\ = 6, v = 1, D = 1. 
K^d presents a maximum only for a > a c ~ 4. Standard cellular conditions (as discussed in the text) correspond to 
green and yellow curves for b/a = 40. c,d, Gain of reactivity due to active transport G±d for d = 1 as a function of 
t 2 (c) and Ti (d) (logarithmic scale). The analytical form obtained in the Methods section (plain lines) is plotted 
against the exact solution (symbols), for the following values of the parameters (arbitrary units): D = l, v = I for 
all curves and a = 10, b = 10 4 (red), a = 10, b — 10 3 (blue), a = 2.5, b = 10 3 (green). Standard cellular conditions (as 
discussed in the text) correspond to blue and red curves. 
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